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Abstract 

We report on theoretical and numerical study of propagation of atomic beams crossing a detuned 
standing-wave laser beam in the geometric optics limit. The interplay between external and internal 
atomic degrees of freedom is used to manipulate the atomic motion along the optical axis by light. By 
adjusting the atom-laser detuning, we demonstrate how to focus, split and scatter atomic beams in a 
real experiment. The novel effect of chaotic scattering of atoms at a regular near-resonant standing 
wave is found numerically and explained qualitatively. Some applications of the effects found are 
discussed. 
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1 Introduction 

Manipulation of atoms by light becomes possible due to the dipole forces which are well described 
by the semiclassical model with quantum description of internal atomic transitions induced by a near 
resonant laser field and classical description of their center-of-mass motion [I]. For the first time, the 
ideas to trap and channel cold atoms with the help of standing laser waves (SLW) have been proposed by 
V. Letokhov and his co-workers [2H3]. The ability of a SLW to deflect, channel and split atomic beams 06] 
has been used for a variety of applications including atom microscopy, interferometry, isotope separation 
and optical lithography [HHS]. Lasers can be used to manipulate atomic trajectories to create atomic 
analogues of such familiar optical phenomena as focusing of light, beam splitting and light scattering. 
It is remarkable that now we are able to reverse the roles of light and matter from their familiar roles. 
The semiclassical description, used in this paper, is similar to the geometric optics limit in conventional 
optics. Atomic trajectories play the role of light rays with the SLW being a light mask. 

In the present paper we intend to demonstrate theoretically and numerically that adjusting in an 
experiment only one parameter, the detuning between the frequencies of a working atomic transition and 
the SLW, one can explore a variety of the regimes of the atom-laser interaction to focus, split and scatter 
atomic beams. Near the atom-field resonance, where the interaction between the internal and external 
atomic degrees of freedom is intense, there is a possibility to create conditions for chaotic scattering of 
atoms [T0HT2] without any additional efforts like a SLW modulation. It becomes possible due to the 
peculiarities of the dipole force in a near-resonant optical lattice |13H16| . 
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2 Focusing, splitting, bunching and scattering of atomic beams 



2.1 Equations of motion 



A beam of two-level atoms in the z direction crosses a SLW laser field with optical axis in the x 
direction. The laser beam has the Gaussian profile exp[— (z — zo) 2 /r 2 ] with r being the e~ 2 radius at 
the laser beam waist. The characteristic length of the atom-field interaction is supposed to be ±1.5r 
because the light intensity at z = ±1.5r is two orders of magnitude smaller than the peak value. The 
longitudinal velocity of atoms, v z , is much larger than their transversal velocity v x and is supposed to 
be constant. Thus, the spatial laser profile may be replaced by the temporal one. The Hamiltonian of a 
two-level atom in the one-dimensional SLW can be written in the frame rotating with the angular laser 
frequency ujf as follows: 

P 2 H 3 ihT 

H = h -(u a ~ Uf)o z - hn exp[-(t - -a t ) 2 /a 2 } (<r_ + &+) cosk f X a — a+6-, (1) 

2m a 2 2 2 

where fr± z are the Pauli operators for the internal atomic degrees of freedom, X a and P are the classical 
atomic position and momentum, T, uj a , and f^o are the decay rate, the atomic transition and maximal 
Rabi frequencies, respectively. The simple wave function for the electronic degree of freedom is (^(t)) = 
a(t)\2) + b(t) 1 1), where a = A + ia and b = B + i/3 are the complex- valued probability amplitudes to find 
the atom in the excited, |2), and ground, |1), states, respectively. 

In the semiclassical approximation, atom with quantized internal dynamics is treated as a point- 
like particle to be described by the Hamilton-Schrodinger equations of motion written for the real and 
imaginary parts of the probability amplitudes 

3 

x = oj r p, ft = — 2exp[— (r — -a T ) 2 / a 2 ](AB + a/3) sinx, 

11 3 

A = -(u r p 2 ~ A)a - -jA - exp[-(r - -cr T ) 2 /a 2 ]/3 cos x, 

113 

« = --j^rV 2 ~ A)^4 - -ja + exp[-(r - -a T ) 2 /al\B cosx, (2) 

1 3 

B = -{uj r p 2 + A)/3 - exp[-(r - -a T ) 2 /a 2 ]a cos x, 

1 3 

P = --{uj r p 2 + A)B + exp[-(r - -a T ) 2 /a 2 ]A cos x, 

where x = kjX a and p = P/Kkf are scaled atomic center-of-mass position and transversal momentum, 
respectively and dot denotes differentiation with respect to the dimensionless time r = £lot. The recoil 
frequency, u) r = Hk 2 /m a fto <C 1, the atom-laser detuning, A = (ujf — u a )/flo, the decay rate 7 = T/Qq, 
and the characteristic interaction time, a T = rQo/v z , are the control parameters. 

Let us introduce instead of the complex-valued probability amplitudes a and b the following real- 
valued variables: 

u = 2Re(ab*), v = -2lm(ab*), z = \a\ 2 - \b\ 2 , (3) 

where u and v are synchronized and quadrature components of the atomic electric dipole moment, 
respectively, and z is the atomic population inversion. In the absence of any losses (7 = 0), Eqs. (|2|) can 
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be cast in the form 



x = oj r p, p=— «exp[— (r a T ) 2 / a 2 } sin x , ii = Av, 

3 3 (4) 

v = — Au + 2exp[— (r — -a T ) 2 /<r 2 ]z cos x, z = — 2exp[— (r — -a T ) 2 /a 2 }v cos x. 

The system @ has two integrals of motion, namely, the total energy 

rr — Wr 2 A , . 

H = —p —ucosx — —2, (5) 

and the length of the Bloch vector, u 2 + v 2 + z 2 = 1, whose conservation follows immediately from 
Eqs. ©. 

Equations @ constitute a nonlinear Hamiltonian autonomous system with two and half degrees of 
freedom which, owing to the two integrals of motion, move on a three-dimensional hypersurface with a 
given energy value H. In general, motion in a three-dimensional phase space in characterized by a positive 
Lyapunov exponent A, a negative exponent equal in magnitude to the positive one, and zero exponent [TT] . 
The maximal Lyapunov exponent characterizes the mean rate of the exponential divergence of initially 
close trajectories and serves as a quantitative measure of dynamical chaos in the system. The values 
of the maximal Lyapunov exponent in dependence on the detuning, the recoil frequency and the initial 
atomic momentum have been computed in Refs. |14|I15|. 

There are different regimes of the center-of-mass motion along the SLW optical axis |12tll5j. In 
dependence on the initial conditions and the values of the control parameters, atoms may oscillate in a 
regular or a chaotic way in wells of the optical potential or move ballistically over its hills with regular 
or chaotic variations of their velocity. Chaotic motion with a positive value of the maximal Lyapunov 
exponent becomes possible in a narrow range of the detuning values, < [A| < 1 |15| . At A = 0, the 
synchronized electric-dipole component, u, becomes a constant. That implies the additional integral of 
motion in the Hamiltonian version (|4|) of Eqs. ([2|) and the regular motion with zero maximal Lyapunov 
exponent. Far off the resonance, at |A| > 1, the motion is regular both in the trapping and ballistic 
modes. 

It is remarkable that there is a specific type of motion, chaotic walking in a deterministic optical 
potential, when atoms can change the direction of motion alternating between flying through the SLW 
and being trapped in its potential wells. We would like to stress that the local instability produces 
chaotic center-of-mass motion in a rigid SLW without any modulation of its parameters. Chaotic walking 
occurs due to the specific behavior of the Bloch-vector component of a moving atom u whose shallow 
oscillations between the SLW nodes are interrupted by sudden jumps with different amplitudes while atom 
crosses each node of the SLW [15]. It looks like a random like shots happened in a fully deterministic 
environment. It follows from the second equation in the set @ that those jumps result in jumps of the 
atomic momentum while crossing a node of the SLW. If the value of the atomic energy is close to the 
separatrix one, the atom after the corresponding jump-like change in p can either overcome the potential 
barrier and leave a potential well or it will be trapped by the well, or it will move as before. The jump-like 
behavior of u is the ultimate reason of chaotic atomic walking along a rigid SW. 

The total atomic energy ([5]) consists of the kinetic one, K = u r p 2 /2, and the potential one, U = 
—ucosx— zA/2. The optical potential changes its depth in course of time. Averaging over fast oscillations 
of the internal atomic variables, we get the averaged potential U = —ucosx — zA/2 that can be used to 
explain why atoms move in such or another way. 
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At small detunings |A| <C 1, the potential is approximately U ~ wcosx. If K(r = 0) > |?7 max | = 1, 
then the atom will move ballistically. This occurs if the initial atomic momentum, pq, satisfies to the 
condition po > \f2/u r . If the initial conditions are chosen to give < K(t = 0) + U(t = 0) < 1, the 
corresponding atoms with < po < ^j2/uj T are expected to move chaotically at the appropriate values 
of A. 

2.2 Focusing and splitting 

In this section we demonstrate how to focus and split atomic beams crossing a Gaussian laser beam 
by varying only one of the control parameters, the atom-field detuning A. Firstly, we perform simulation 
with a negligible probability of spontaneous emission and solve the Hamiltonian equations of motion §4§ 
at comparatively large value of the detuning, A = 1. To be concrete we take as an example calcium 
atoms with the working intercombination transition 4 1 S'o — 4 3 Pi at A a = 657.5 nm, the recoil frequency 
f rec ~ 10 KHz, and the lifetime of the excited state T sp = 0.4 ms. Taking the maximal Rabi frequency 
to be Q,q/2tt = 2 • 10 7 Hz, the radius of the laser beam r = 0.3 cm, and the mean longitudinal velocity 
v z = 10 3 m/s, the interaction time is estimated to be 0.9 ms, longer than the atomic lifetime. The 
normalized recoil frequency is oj r = 47Tf r ec/^o = 10 -3 and the normalized characteristic time is a T = 400. 

Trajectories for 50 calcium atoms to be prepared in the ground states (uq = vq = 0, zq = —1) with 
the same initial momentum, pq = 10, and initial positions in the range — 7r/10 < x < n/W are shown in 
Fig.[TJ In units of the optical wavelength, X = x/2ir, this range is —0.05 < X < 0.05. The focusing occurs 
at those moments of time when the average transverse momentum in the atomic beam is approximately 
equal to zero. If one turns off the laser at one of these moments, it becomes possible to reduce the beam 
width practically in ten times. The reason of focusing is simple. It is well known [Jj that at positive blue 
detunings atoms are attracted to the nodes of the SLW where the minima of the optical potential are 
situated at A > 0. The first node the atoms reach at po > is situated at X = 1/4. The initial kinetic 
atomic energy, Kq = 0.05, is not enough to overcome the potential barrier whose depth can be estimated 
to be ~ 0.35 because the simulation gives u ~ and z ~ —0.7. So, all the atoms in the beam oscillate in 
the first potential well in the x-direction around the first node. The initial width of the beam, 5Xo = 0.1, 
is gradually reduced because in course of time the atoms with initial negative positions catch up with the 
ones with positive Xq near the first turning point where the average beam momentum is close to zero. 
The time interval of the atomic interaction with the SLW field is estimated to be 3a T = 1200. So, the 
atoms leave the potential well after that time and move freely (see Fig. [1]) . 

In order to take into account spontaneous emission we use the standard stochastic wave-function 
technique [18H20j for solving Eqs. (|2|). The integration time is divided into a large number of small 
time intervals St. At the end of the first interval, r = r%, the probability of spontaneous emission, 
si = j5t\o, Ti | 2 /(|a Tl | 2 + l^nl 2 ); is computed and compared with a random number, e, from the interval 
[0, 1]. If s\ < ei, then one prolongs the integration but renormalizes the state vector in the end of the 
first interval at r = : a T + = a Tl /^J\a Tl | 2 + \b Tl \ 2 and b T + = b Tl / \f\a Tl | 2 + \b Tl | 2 . If si > £i, then the 
atom emits a spontaneous photon and jumps to the ground state at r = r\ with A T1 = a T1 = f3 T1 = 0, 
B Tl =1. Its momentum in the x direction changes for a random number from the interval [0, 1] due to 
the photon recoil effect, and the next time step commences. 

We simulate lithium atoms with the relevant transition 25 , ]/ 2 — 2P 3 / 2 , the corresponding wavelength 
A a = 670.7 nm, recoil frequency f rec = 63 KHz, and the decay time T sp = 2.73- 10 -8 s. With the maximal 
Rabi frequency f2o/27r ~ 126 MHz and the radius of the laser beam r = 0.05 cm one gets io r = 10~ 3 , 
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Figure 1: (a) Focusing the atomic beam with a long lifetime of the excited state, (b) The effect of 
spontaneous emission on the focusing. The detuning is A = 1 in both the cases. The atomic position X 
is in units of the optical wavelength. 




a T = 400, and 7 = 0.05. Trajectories for 50 spontaneously emitting atoms under the same conditions as 
in Fig. are shown in Fig. [lb. As expected, spontaneous emission destroys in part the effect of focusing. 
However, the atoms move more or less coherently because spontaneous emission events are comparatively 
rare at A = 1. 

The other effect, we would like to demonstrate with atomic beams crossing the SLW, is a splitting of 
the beam. To do this one needs to choose such the value of the detuning in order that some atoms in 
the beam would be trapped in the first well of the optical potential but another ones could overcome the 
barrier and leave that well. It is possible to split atomic beams as at positive and negative values of the 
detuning. As an example, we demonstrate in Fig. [2] the effect of splitting at A = —1 for atoms without 
and with spontaneous emission. It is seen that spontaneous emission changes slightly the effect because 
a few atoms may leave the potential wells due to random recoils. 
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Figure 3: The effect of bunching of 10 5 calcium atoms around (a) the SLW nodes (blue detuning, A = 1) 
and (b) the antinodes (red detuning, A = —0.2). The plot of atomic density n = N(X)/Nq at the fixed 
moment of time r = 1000. 



3 Bunching and chaotic scattering of atoms 



The ability of blue and red detuned lasers to attract atoms to the nodes and antinodes of the SLW, 
respectively, can be used to create periodic structures composed of atoms deposited on substrates in 
the process of optical nanolithography [THS]- To simulate a real experiment we consider a beam with 
Nq = 10 5 calcium atoms with the initial Gaussian distribution (with the rms a x = a v = 2 and the average 
values xo = and po = 10) and compute their distribution against the SLW at a fixed moment of time 
r = 1000. The bunching of atoms at the SLW nodes at A = 1 (blue detuning) is shown in Fig. [3k where 
the atomic density, n = N(X)/Nq, is plotted along the optical axis X at r = 1000. The same effect, but 
with the atoms bunching around the SLW antinodes (red detuning, A = —0.2), is shown in Fig. [3]3. In 
both the cases we get a periodic atomic relief with the period A//2 the width of which is restricted by 
the time the atoms interact with the Gaussian laser beam. 

The problem we consider resembles the scattering process with particles entering an interaction re- 
gion along completely regular trajectories and leaving it along asymptotically regular trajectories. It is 
known from many studies in celestial mechanics, fluid dynamics and other disciplines that under certain 
conditions the motion inside the interaction region may have features that are typical for dynamical 
chaos, (homoclinic and heteroclinic tangles, fractals, strange invariant sets, positive finite-time Lyapunov 
exponents, etc.) although the particle's trajectories are not chaotic in a rigorous sense because chaos 
is strictly defined as an irregular motion over infinite time. It has been found |21H24j that transient 
Hamiltonian chaos in the interaction region occurs due to existence of, at least, one nonattractive chaotic 
invariant set consisting of an infinite number of localized unstable periodic orbits and aperiodic orbits. 
This set possesses stable and unstable manifolds extending in the phase space into the regions of regular 
motion. The particles with the initial positions close to the stable manifold follow the chaotic-set tra- 
jectories for a comparatively long time, then deviate from them, and leave the interaction region along 
the unstable manifold. It is the common mechanism of chaotic scattering that in our problem causes the 
chaotic walking of atoms along the SLW. 

In Fig. [3k we show the atomic position distribution at r = 1000 in the regime of the chaotic scattering 
at A = 0.2 with 10 5 calcium atoms. This plot should be compared with Fig. [3k where the atomic position 



6 



n*io 

6 



a) 









mo' 2 



-5 -4 -3 -2 -1 1 2 3 4 5 




30 p 



Figure 4: The distributions of 10 5 calcium atoms at r = 1000 in (a) the real and (b) momentum space 
under the conditions of chaotic scattering at A = 0.2. 





Figure 5: Comparison of the distributions of 10 5 calcium atoms at r = 1000 over the phase plane in the 
regimes of (a) chaotic (A = 0.2) and (b) regular scattering (A = —0.2). 
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distribution is shown for regularly scattered atoms at A = 1. First of all, the distribution of chaotically 
scattered atoms has a prominent pedestal and is much broader. Moreover, it has no such a periodic 
structure as shown in Fig. [3^. Only the peaks around the first two SLW nodes are prominent in Fig. 0k. 
The atomic position distribution in the momentum space in Fig. is much broader than the one for 
regularly scattered atoms at A = 1 (not shown). Thus, we predict that under the conditions of chaotic 
scattering there should appear less contrast and more broadened atomic reliefs as compared to the case 
of regular scattering because a large number of atoms are expected to be deposited between the nodes 
as a result of chaotic walking along the SLW axis. The effect is expected to be more prominent under 
the coherent evolution but it seems to be observable with spontaneously emitting atoms as well. The 
difference between chaotic and regular scattering of atoms at a rigid SLW is especially prominent on the 
corresponding phase space portraits shown in Fig. [5] where positions and momenta of 10 5 calcium atoms 
are plotted at the fixed time moment. 



4 Conclusion 



We have simulated some geometric optics effects with atomic beams crossing a SLW in the limit of long 
relaxation time and with spontaneous emission taken into account. Trajectories of spontaneously emitting 
atoms have been simulated with the help of the standard stochastic wave-function technique [18H20j . It 
has been shown that by adjusting the detuning it is possible to focus, split and scatter atoms. The 
effects have been explained by a coupling between external and internal atomic degrees of freedom. The 
depth of the optical potential depends on the sign and value of the detuning. Varying A, one can create 
conditions for focusing, splitting and bunching the atoms. It is remarkable that near the atom-field 
resonance we have found the new type of atomic diffraction at a SLW without any modulation of its 
parameters that can be observed in real experiments. That would be the prove of existence of the novel 
type of atomic motion, chaotic walking in a deterministic environment. The effects found could be used 
in optical nanolithography to fabricate complex atomic structures on substrates. 

We predict that experiments on the scattering of atomic beams at a SLW can directly image chaotic 
walking of atoms along the SLW. In a real experiment the final spatial distribution can be recorded 
via fluorescence or absorption imaging on a CCD, commonly used methods in atom optics experiments 
yielding information on the number of atoms and the cloud's spatial size. The other possibility is a 
nanofabrication where the atoms after the interaction with the SW are deposited on a silicon substrate 
in a high vacuum chamber. In this case the spatial distribution can be analyzed with an atomic force 
microscope. As to the momentum distribution, it can be measured, for example, by a time-of-flight 
technique [25J. The modern tools of atom optics enable to create narrow initial atomic distributions in 
position and momentum, reduce coupling to the environment and technical noise, create one-dimensional 
optical potentials, and to measure spatial and momentum distributions with high sensitivity and accuracy. 



Acknowledgments 

This work was supported by the Russian Foundation for Basic Research (projects nos. 09-02-00358 
and 09-02-01258), by the Integration grant from the Far-Eastern and Siberian branches of the Russian 
Academy of Sciences, and by the Program "Fundamental Problems of Nonlinear Dynamics" . 



8 



References 

[1] A.P. Kazantsev,G.A. Ryabenko,G.I. Surdutovich, V.P. Yakovlev, Phys. Rep., 129, 75 (1985). 

[2] V.S. Letokhov, JETP Lett, 7, 272 (1968) [Pis'ma ZhETF, 7, 348 (1968)]. 

[3] V.I. Balykin, V.S. Letokhov, Opt Comm., 64, 151 (1987). 

[4] V. Letokhov, Laser control of atoms and molecules (Oxford University Press, New York, 2007). 

[5] E. Arimondo, A. Bambini, S. Stenholm, Phys. Rev., 24, 898 (1981). 

[6] C.S. Adams, M. Sigel, J. Mlynek, Phys. Rep., 240, 143 (1994). 

[7] T. Sleator, T. Pfau, V. Balykin, O. Carnal et al, Phys. Rev. Lett, 68, 1996 (1992). 

[8] G. Timp, R.E. Behringer, D.M. Tennant, J.E. Cunningham, Phys. Rev. Lett, 69, 1636 (1992). 

[9] J.J. McClelland, R.E. Scholten, E.C. Palm, R.J. Celotta, Science, 262, 877 (1993). 

[10] S. V. Prants, L.E. Kon'kov, JETP Letters, 73, 1801 (2001) [Pis'ma ZhETF, 73, 200 (2001)]. 

[11] S. V. Prants, M. Edelman, G. M. Zaslavsky, Phys. Rev. E, 66, art. 046222 (2002). 

[12] S.V. Prants, V.Yu. Sirotkin, Phys. Rev. A, 64, 033412 (2001). 

[13] V. Yu. Argonov, S. V. Prants, JETP, 96, 832 (2003) [ZhETF, 123, 946 (2003)]. 

[14] V. Yu. Argonov, S. V. Prants, J. Russ. Laser Res., 27, 360 (2006) 

[15] V. Yu. Argonov, S. V. Prants, Phys. Rev. A, 75, art. 063428 (2007). 

[16] V. Yu. Argonov, S. V. Prants, Phys. Rev. A, 78, art. 043413 (2008). 

[17] L.E. Kon'kov, S. V. Prants, JETP Letters, 65, 833 (1997) [Pis'ma ZhETF, 65, 801 (1997)]. 

[18] H. J. Carmichael, An open systems approach to quantum optics (Berlin, Springer, 1993). 

[19] J. Dalibard, Y. Castin, K. Molmer, Phys. Rev. Lett, 68, 580 (1992). 

[20] R. Dum, P. Zoller, H. Ritsch, Phys. Rev. A, 45, 4879 (1992). 

[21] P. Gaspard, Chaos, Scattering and Statistical Mechanics, Cambridge University Press, Cambridge 
(1998). 

[22] K.A. Mitchell, J.P Handley, B. Tighe, J.B. Delos, S.K. Knudson, Chaos, 13, 880 (2003). 

[23] M. Budyansky, M. Uleysky, S. Prants, Physica D, 195, 369 (2004). 

[24] M.V. Budyansky, MYu. Uleysky, S.V. Prants, JETP, 99, 1018 (2004) [ZhETF, 126, 1167 (2004)]. 

[25] M. G. Raizen, Adv. At. Mol. Opt. Phys., 41, 43 (1999). 



9 



